Theory of Myelin Coiling 
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Abstract. A new model is proposed to explain coiling of myelins composed of fluid bilayers. This 
model allows the constituent bilayer cylinders of a myelin to be non-coaxial and the bilayer lateral 
tension to vary from bilayer to bilayer. The calculations show that a myelin would bend or coil to lower 
its free energy when the bilayer lateral tension is sufficiently large. From a mechanical point of view, 
the proposed coiling mechanism is analogous to the classical Euler buckling of a thin elastic rod under 
axial compression. The analysis of a simple two-bilayer case suggests that a bilayer lateral tension of 
about 1 dyne/cm can easily induce coiling of myelins of typical lipid bilayers. This model signifies the 
importance of bilayer lateral tension in determining the morphology of myelinic structures. 

PACS. 87.16.Dg Membranes, bilayers, and vesicles - 82.70.Uv Surfactants, micellar solutions, vesicles, 
lamellae, amphiphilic systems, etc. - 82.70.-y Disperse systems; complex fluids 



1 Introduction 

Surfactant molecules such as polar lipids self-assemble into 
fluid bilayers when dissolved in water at sufficiently high 
temperatures PE]. If the surfactant concentration is large 
enough, bilayers stack to form multilayer structures called 
lamellae or multilamellae {La phase) 0. In some ex- 
periments ^ El El El El El) bilayers curve collectively to 
form nested cylindrical multilamellae known as myelin fig- 
ures or simply, myelins (Fig. Under certain condi- 
tions, myelins may bend to form coils or double helices 
IHllTnillllinilTnillll. Because myelins have potential ap- 
plications such as controlled drug delivery detailed 
knowledge of their structure and behavior is desired. Re- 
cently, a geometrical model was proposed to account for 
formation of myelins ■ This model suggests that if the 
inter-bilayer repulsion is large enough, a flat multilamella 
is unstable against myelin formation, i.e., it would trans- 
form into a myelin or myelins in order to lower its energy. 
In this paper we will focus on coiling of myelins composed 
of fluid bilayers and present a theory to explain why they 
coil. 

Several theories of myelin coiling have been proposed. 
Lin et al. reported that the presence of enough Ca^+ 
ions led to coiling of the myelins of binary mixtures of 
cardiolipin and phosphatidylcholines Mishima and 

his co-workers studied the double helix formation of the 
myelins composed of egg-yolk phosphatidylcholine bilay- 
ers 22 ■ The authors of these two works attributed the 




Figure 1: Sketch of myelin structure. Myelins (M) are 
multilamellar tubes composed of nested bilayer cylinders. 
Oval inset shows arrangements of surfactant molecules (L) 
in bilayers (B) ^E], where the mid-surface (S) of a bilayer 
separates the two opposing monolayers. The symbol d de- 
notes the inter-bilayer spacing. In the upper drawing (M), 
only the bilayer mid-surfaces are shown. Typical value of 
bilayer thickness or inter-bilayer spacing is about several 
nanometers (l, ,2, and that of myelin diameter ranges from 
a few micrometers to about 50 micrometers |S1 El|3 El El . 
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observed coiling to surface adhesion. Frette et al. investi- 
gated coiling of the myelins of phospholipid bilayers in the 
presence of anchored polymers They suggested that 
coiling was due to a coupling between local bilayer curva- 
ture and polymer concentration. In addition to the above 
theories, Santangelo and Pincus proposed a myelin model 
and showed that myelins are unstable to coiling when the 
bilayer spontaneous curvature is nonzero or when the equi- 
librium inter-bilayer spacing is decreased |17| . Here we de- 
velop a new model. Our model allows the bilayer cylinders 
of a myelin to be non-coaxial and the bilayer tension to 
vary from bilayer to bilayer ^j. By taking account 
of these two degrees of freedom, we will show that a suf- 
ficiently large bilayer tension may cause a myelin to coil. 
The coiling mechanism proposed here is similar to the clas- 
sical Euler buckling of a thin elastic rod under axial load 

This paper is organized as follows. In Section |21 we 
describe the experimental observations that lead to our 
model. The proposed coiling mechanism is explained with 
a two-dimensional example in Section 13 The definition of 
our myelin coiling model is given in Sectional We study 
coiling of a myelin composed of multiple bilayers in Section 
|S1 Coiling of a two-bilayer myelin is investigated closely in 
Section 1^1 In Section |7| we discuss important features and 
implications of our model. Finally, Section|Slconcludes our 
work. 



Contact Experiment and Myelin 
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In order to provide the motivation for this work, we will 
examine closely how myelins are created as well as their 
internal structure. We will show that the lateral tension 
of bilayers |18[ [T^ should be taken into consideration and 
moreover, the constituent bilayer cylinders of a myelin need 
not be concentric. 

Myelins are multilamellar tubes composed of nested 
cylindrical fluid bilayers (Fig. They are usually pro- 
duced in contact experiments |S] : water is brought 
into contact with a lump of concentrated surfactant at 
temperatures high enough so that the bilayers are in fluid 
state. Upon contact myelins form at the interface and grow 
into water. The growth of myelins is roughly diffusion-like 
initially and slows down significantly after a few minutes 
I14L 1^ . While myelins are still growing, coiling of sin- 
gle myelins or double helix formation of two intertwisting 
myelins are often observed JT]. Experiments show that 
as myelins grow, water enters myelins mainly via the root 
region (Fig. [7||H|. Thus, the influx of material, surfac- 
tant and water, that supports myelin growth comes from 
the back of myelins, i.e., the root region. Because the root 
region contains defects of sizes smaller than typical myelin 
diameters [^IIISIj the material flux density for a growing 



Figure 2: Sketch of myelin growth and material flux in a 
contact experiment, (a): Myelins (M) grow from the inter- 
face (IF) of water (W) and bulk surfactant (S). The wa- 
ter flux (WF) and surfactant flux (SF) supporting myelin 
growth enter myelins via the root region (RR). Solid and 
dash arrows indicate the directions of the water and sur- 
factant fluxes, respectively [7||S|. (b): A single myehn (M) 
is naively divided into a core (I) and a shell (O). As repre- 
sented by the lengths of the respective arrows, the material 
flux densities for these two parts can be different. The bi- 
layers of the shell would be in tension if the material flux 
density for the core is larger than that for the shell. 



2 



myelin is expected to be nonuniform^ . Figurel^b provides a 
simple and idealized example of nonuniform material influx 
density. In this example the myelin is arbitrarily divided 
into two parts, a core and a shell. The defect structure in 
the root region is such that the material flux density for 
the shell is less than that for the core, which implies that 
the shell is being stretched axially and the core compressed 
also axially. Therefore the bilayer tension is expected to 
be nonuniform and furthermore, the bilayers in the shell 
should be under tension. We note that the fluidity of a bi- 
layer leads to a uniform tension throughout the bilayer. In 
contrast, the notion of "nonuniform tension" used in this 
work always means the tension may vary from bilayer to 
bilayer. 

Intuitively the bilayer tension in a myelin composed of 
fluid bilayers should have some degree of nonuniformity, 
regardless of how the myelin is prepared. Like an air bubble 
in water, where the interfacial tension is balanced by a 
pressure difference across the air-water interface (Laplace's 
formula) , at least the outermost bilayer of a myelin should 
be under tension in order to cancel out the repulsion from 
the neighboring bilayer and thus maintain its cylindrical 
shape of the body and spherical shape of the end cap(s). 
If one or some of the bilayers are under tension, the force 
balance for the myelin body or the end cap(s) implies that 
the bilayer tension may be nonuniform, i.e., it may change 
from bilayer to bilayer. 

From the above arguments we expect that the bilayer lat- 
eral tension can be variable not only from bilayer to bilayer 
but also from myelin to myelin. A general myelin model 
should take this degree of freedom into account. Myelins 
of uniform bilayer tension are just special cases of such a 
model. In addition, we will not presume the bilayer cylin- 
ders of a myelin to be concentric or coaxial since there is no 
reason to do so, especially for coiled myelins. Hence, the 
eccentricities of the bilayer cylinders appear in our myelin 
model naturally. We will show how these two variables, 
the bilayer lateral tension and cylinder eccentricity, can 
give rise to myelin coiling. 

3 Bending of a Two-dimensional 
Tube 

This section gives a detailed account of how an artifi- 
cial two-dimensional tube can be unstable against bend- 
ing, i.e., bending may reduce the tube's free energy. Al- 
though myelin coiling occurs in three dimensions, this two- 
dimensional case, nevertheless, captures all the essential 
features of the myelin coiling model proposed in this work. 
Therefore, not only does this simple example offer a phys- 
ical insight into why myelins coil, it also provides a frame- 

^The nonuniformity of the material flux density depends on the 
defect structure in the root region, which would evolve over time 
since defects are not energetically favorable. 



work for modeling myelin coiling later in this paper. In the 
following we will first define the two-dimensional tube and 
specify how it deforms. Then we will show that the tube 
may indeed curve to lower its energy. The two-step defor- 
mation developed here makes it easy to determine changes 
in the energies. 

3.1 Geometry and Two-step Deformation 

Figure |31 describes the geometry of the two-dimensional 
tube and explains how it deforms. The tube comprises a 
core, fc, and a perimeter, abde, which are two-dimensional 
analogues of the inner and outer cylinders of the two- 
bilayer myelin studied in Section respectively. All the 
state parameters, such as dr, di and 9, are defined in the 
figure caption. The eccentricities of the bilayer cylinders 
in a myelin are simulated in this two-dimensional case by 
allowing the spacings di and di- between the core and the 
perimeter to vary. The perimeter abde of length 

L2D =Li + L^ + L, (1) 

is extensible and under uniform tension. The core fc exerts 
a force on the end cap bd to counterbalance the tension of 
the perimeter and prevent the tube from collapsing. For 
simplicity we assume the core to be of zero thickness and 
subject to bending only, i.e., its length is constant. 
One may think of the core as a spring of infinite spring 
constant. Allowing to vary amounts to reducing the 
spring constant to a finite value and thus adding one more 
degree of freedom to the system, which would not invali- 
date the conclusions of this section. 

The tube area A2B is postulated to be constant (Fig. |3J|: 

A2D = 2rfo-^o + ^dl ■ c 

= ^-[{R+d,f-iR-d,r] 

+ |(^)%, (2) 

where the angle 9 can be larger than 2tt, i.e., the tube 
may coil, and c = or 1. The area conservation in two di- 
mensions is analogous to the volume conservation in three 
dimensions (Sect. I4.3|l . In order to examine the effect 
of the end cap bd, we deliberately multiply every end cap 
contribution by the parameter c. We will demonstrate that 
when the tube's aspect ratio Lo/2do is much greater than 
one, all the end cap contributions are negligible, i.e., we 
can set c = 0. 

The deformation of this two-dimensional tube takes 
place in two steps (Fig. |3Jl: the tube is transformed from 
the initial state (i) into the intermediate state (ii) and 
then into the final state (iii). The initial tube is straight 
and its shape is symmetric with respect to the line fc. 
From state (i) to state (ii) the tube curves with fixed di, 
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Figure 3: Two-step deformation of a two-dimensional 
tube. The tube is initially in state (i) and transforms into 
the intermediate state (ii) and then into the final state (iii) . 
During the deformation the tube area is conserved. The 
tube consists of a core, fc, and a perimeter, abde. The 
core fc is of zero thickness and subject to bending only, 
i.e., its length is fixed. The perimeter abde is composed 
of two lines (arcs), ab and de, and a half circle, bd. As 
indicated in the figure, the symbols L\, Lc, Lr, L^^, di and 
dr denote the lengths of the lines (arcs) ab, bd, de, fc, af 
and fe, respectively. Initial state (i): The tube is straight 
and its shape is symmetric with respect to the line fc, i.e., 
di = dr = do, L\ = Lj- = = Lq and thus Lc = Trdo- 
Intermediate state (ii): The tube curves with d\ = d-^ = do 
and — R9 = Lq. The point O is the center of bending 
curvature, R the radius of curvature of the arc fc, and 9 
the angle subtended by fc. The lengths Ly = {R — do)0, 
Lr — {R + do)d, and Lc = Trdo- Final state (iii): The 
spacing d-^ decreases from do with fixed R and 9. In state 
(iii), di > do > dr, im = R9 = Lq, L\ = {R - d\)9, 
Lr — {R + dr)9, and Lc = n{di + dr)/2. Because of area 
conservation |5J), the width {d\ + dr) of state (iii) is greater 
than the initial value 2do. 



dr and L^. Consequently, the perimeter length L2D iQl 
is also unchanged. In the final state (iii) the spacing dr 
decreases from its original value do while Lm is still un- 
changed. Equation (|2Jl gives the change in the perimeter 
length L2D from state (i) or (ii) to state (iii) 
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R ' 2Lo 

where Ad = di — dr. It is obvious that SL2D < when 
Ad > and do ^ Lq. In state (iii) Ad has to be greater 
than zero because of the area conservation (|2J). We now 
investigate how the tube's free energy changes with the 
parameters R and Ad- 

3.2 Free Energy and Instability 

The total free energy of the tube is given by 



^2D 



-^213 



2D J 



(4) 

where -F^p is the bending energy and the elastic energy 
associated with the tension in the perimeter. The interac- 
tion between the core and the perimeter is neglected for 
the sake of simplicity as well as clarity. We assume that 
the bending energy F]'^ takes the form of 

K2D 



^2D 



i?2 



■Lq + aoc, 



(5) 



where K2d is the bending stiffness of the tube, R the radius 
of bending curvature, and Lq = R9 the tube length (Fig. 
131). The last term oqc represents the contribution of the 
end cap bd. The tension E in the perimeter is defined as 
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(6) 



where L2D is the perimeter length Now we can deter- 
mine the difference in the total free energy F2B between 
the initial state (i) and the final state (iii). 

When the tube is deformed from state (i) to state (ii), 
only the bending energy F]}-, © changes because di, dr, 
L,„ and thus L21) arc unchanged (Fig. |2J|. Therefore 

■2D[(ii)]-^^2'D[(i)] ^'"^ 



f:? 



i?2 



■Lo. 



(7) 



The change in F^^^ from state (ii), where Ad = 0, to state 
(iii), where Ad > 0, is evidently given by 



J^2D[(iii)]-^2D[(n)] = 



h{R) 
R 



LoAd + b2c + 0{Al). (8) 



The first two terms on the right-hand side of the above 
equation represent the lowest-order correction of the bend- 
ing energy due to a change in Ad. The term b2C, also a 
function of Ad, is the end cap contribution. Combining Q 
and IHl) leads to 
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^2dI 
K2D 



(iii)] 



-'^2D 

bliR) 



(i)] 



-^Lo + -^LoA^ + b2C + OiAi) (9) 
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For cases with finite R, smah Ad, and large aspect ratio 
{Lq ^ 2do), equation ||HJ) is much less than equation iQ 
and thus can be neglected. Equations (0) and © give the 
change in the elastic energy -F2D- 

^2D[(iii)]-^2D[(i)] 
1 OS 



E • 6L2D 



R 



-1 



2 9i2D 

ttcIqc 
~2L^ 



{SL2Br + ■ 
0{Al). 



(10) 



Equations © and (|10|l lead to the change in Q 
from state (i) to state (iii): 
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where &iAd is much less than k,2d/R and thus neglected 
(see the text below 0). Ahhough 6F^j^ > 0, SF^j^ and 
hence dF2B can be negative. Consequently, if the perimeter 
tension S is large enough, a tube in state (i) would curve 
to lower its free energy. The equilibrium value of {R, Ad) 
can be determined if the 0{A'^) term in is known. 
In addition, all the above equations involving the end cap 
show that as long as the tube is sufficiently long, the end 
cap contributions are negligible, i.e., c can be set to zero. 

3.3 Summary 

We have shown that a two-dimensional tube of a given 
area is unstable against bending when the perimeter ten- 
sion is sufficiently large. Owing to the area conservation, 
the perimeter length decreases when the tube undergoes 
the two-step deformation (Fig. The decrease in the 
perimeter length reduces the elastic energy F|j-, © . There- 
fore, the tube can bend to lower its energy if the reduction 
in F2Y) is larger than the increase in the bending energy 
^2D- ^Iso have shown that the end cap energy and its 
changes are negligible when the tube has a large aspect 
ratio, i.e., Lo/2do » 1. 

This two-dimensional case will serve as a paradigm for 
modeling myelin coiling in the following sections. In ad- 
dition to the bilayer bending rigidity and the bilayer lat- 
eral tension, the inter-bilayer interaction will also be taken 
into consideration in our myelin model. We will employ 
the two-step deformation (Fig. to determine the energy 
changes associated with myelin coiling. 

4 Model Definition 

The myelin model proposed in this section takes the bi- 
layer tension and cylinder eccentricity into account (Sect. 
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Figure 4: Model of an iV-bilayer myelin. The myelin is 
composed of N nested bilayer cylinders. These cylinders 
may not be coaxial and moreover, they can bend coop- 
eratively to yield myelin coiling (Figs. |S1 and IHJ). The 
symbol r„ denotes the radius of the rt-th cylinder, where 
1 < n < A^. Since these bilayer cylinders need not be 
concentric, we label their diameters 2r„ instead of their 
radii r„. The length L„ of the n-th bilayer cylinder is de- 
fined as that of its axis (Fig. [7|l, and the length Lm of a 
myelin as that of the innermost cylinder, i.e., Lm = Li. 
For a straight myelin, L„ and Lm are equal for all n, 
whereas for a curved myelin, they may not (Fig. (HJ. In 
this work we only consider myelins of large aspect ratios, 
i.e., Lm/^tm ^ 1. Therefore the end caps (EC) of myelins 
can be neglected. 

121) • We postulate that the growth of myelins is sufficiently 
slow so that the myelins are in quasi-equilibrium and there- 
fore their energies are well defined In the following 
we will first specify the geometry and energy of a myelin 
and then describe how myelin coiling is modeled in this 
work. Like the two-dimensional case of the previous sec- 
tion, the proposed coiling mechanism is analogous to the 
Euler buckling of an elastic rod subject to axial compres- 
sion 1201 • 

4.1 Myelin Geometry 

An A-bilayer myelin is a multilamellar tube composed of a 
nested stack of A^ cylindrical fluid bilayers (Fig. . These 
bilayer cylinders need not be coaxial (Fig. and may 
curve collectively to give rise to myelin coiling (Fig. inj. In 
addition, we only consider myelins of large aspect ratios so 
that the end caps are negligible. The symbol r„ denotes the 
radius of the n-th bilayer cylinder, with n = 1 and n = N 
corresponding to the innermost and outermost cylinders, 
respectively (Fig. 0J. We define the length L„ of the n-th 
bilayer cylinder as that of its axis and the length Lm of 
a myelin as that of the innermost cylinder, i.e., Lm = Li 
(Figs. 2] and |7|). With this myelin model we can now 
describe the geometry of a coiled myelin. 

A coiled myelin is composed of nested helical bilayer 
cylinders (Fig. The following parametric expression 
for a helix describes the axis of the n-th bilayer cylinder: 

R„((/)) = (p„cos(/), p„sin(/), P0/27r), (12) 

where 0, p„ and P are defined in Figure El The radius of 
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Figure 6: Coiling of a myelin. For clarity we only show 
the n-th bilayer cylinder (Fig. 0J|. The axis of this coiled 
cylinder forms a helix described by the parametric equation 
Rn,((/)) — [pn cos (j), pn sin (j), P0/27r), where pn is the radius 
of the circular projection of the helix on the X-Y plane and 
P the pitch of the helix. 



Figure 5: Cross-sectional view of two adjacent bilayer 
cylinders in a myelin. We describe the n-th and [n + 1)- 
th cylinders with the vectors S„ and S„+i, and their 
axes with 0„ and 0„+i, respectively. When two nested 
cylinders are not concentric, the spacing between them 
is clearly not uniform (Fig. (HJ. In this work we define 
the spacing d„ between the n-th and (n + l)-th cylinders 
as d„(S„) = min{|S„+i — S„|}. To calculate (i„(S„), we 
select the coordinate system with 0„ being the coordi- 
nate origin and S„ — (r„ cosV-'n, sin V'n)- The eccen- 
tricity vector E„_|_i = {Sn+i,0) represents the positional 
difference between On+i and 0„, as shown in the oval 
inset. Given Sn+i, the inter-bilayer spacing (i„(S„) = 
r„+i — (r^j -I- Sj^_^_i — 2r„(5„+i cos V'n)^^^- For convenience 
we set i5i = 0. 













\ ^ \ 
V' ' 





Figure 7: Bending of the n-th bilayer cylinder. The point 
O is the center of bending curvature. The cylinder length 
Ln = RnQ, where i?„ is the radius of bending curvature 
H13() and 9 the angle subtended by the cylinder. This 
cylinder can be viewed as a section of a torus of radii 
r„ and i?„. Straight cylinders correspond to cases with 
Rn — > oo, 6* — > 0, and finite i?„6'. The surface area element 
dAn — rnd^/jn{Rn + costpn)d9. The two principal curva- 
tures of the surface are l/r„ and cosipn/ {Rn + fn cosipn)- 
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curvature of the helix R„ is equal to 

p2 



Rn — Pn 1 



(13) 



The volume Vn, area A„ and length L„ of the n-th bilayer 
cylinder are given by 



Vn = 7rr„L„, 



27r/9„ 



p2 



47r2p2 



1/2 



(14) 
(15) 

(16) 



where r„ is the radius of the cylinder and the number of 
turns of the coiling (Fig. Strictly speaking, An is the 
projected area of the n-th bilayer J^l- The true bilayer 
area is larger than An because of thermal undulations of 
the bilayer. 

4.2 Free Energy 

The total energy of a myelin includes the curvature en- 
ergy F'^ of bilayers, the elastic energy F° associated with 
bilayer lateral tension, and the inter-bilayer interaction en- 
ergy F': 

F = F"+F° + F'. (17) 

The curvature energy E'^ of a membrane of area A takes 
the form of 



^ - 2 



dA{K^ +K2- Ko) 



(18) 



where k is the bending stiffness, and Ki and K2 are the 
principal curvatures |25|. The spontaneous curvature Kq 
is taken to be zero for simplicity. (See the Discussion sec- 
tion.) The Gaussian curvature term (~ J dAKiK2) is ne- 
glected because the myelin coiling investigated here does 
not involve any topological changes. Thus, there can be 
no change in this energy according to the Gauss-Bonnet 
theorem [21]. Because of the fluid nature of bilayers, a 
myelin can not sustain torsional stress in its axial direction. 
Therefore the curvature energy change of a myelin due to 
coiling depends only on the bending of its constituent cylin- 
ders (Figs. El and 13). Using lO and (0, with the area 
element dAn = rndipn{Rn + cos'0n)d^ and the two prin- 
cipal curvatures 1 / r„ and cos ipn/ (Rn + fn cos ipn) (Fig. , 
the curvature energy E!^ of the n-th bilayer cylinder in a 
coiled myelin is given by 



E?, = 



TxnLn 



{rn/Rn) 



2RI 



, 3 rl 
1 -I — 

iRl 



(19) 



where Z/„ = RnO is the cylinder length. The above equa- 
tion suggests that the curvature energy of a curved bi- 
layer cylinder can be written as a sum of the energy 



(~ l/j'n) associated with local bending of the bilayer and 
that (~ l/Rn) resulting from bending of the cyhnder fl^ . 
Equation lfll?|l leads to the curvature energy of a coiled 
A^-bilayer myelin: 



N 

E 

n=l 
N 

E 

n=l 



E'„ 



N ^ 



(20) 



where the bending modulus k„ of the n-th bilayer cylinder 
is defined as 



TTKrn 



1 



3< 



-I- 



(21) 



The lateral tension of a bilayer arises from stretching 
the bilayer area from the relaxed state |19[ I25L ESj • For a 
fluid bilayer, the tension is uniform and isotropic. Let F^ 
denote the elastic energy giving rise to the tension cr„ in 
the n-th bilayer cylinder, i.e.. 



dAn 



(22) 



where An is the area of the cylinder H15(l . The elastic 
energy F° of an A^-bilayer myelin is then given by 



N 



(23) 



ra=l 



The energy describes the net result of all kinds of 
inter-bilayer interactions such as van der Waals attraction, 
electrostatic repulsion [2] , hydration pressure 1^, and the 
Helfrich repulsion [2H|- For an A^-bilayer myelin, F' can be 
conveniently expressed as 



N-l 

F' = V Fl 



(24) 



where F^ represents the energy of interaction between the 
n-th and (n -I- l)-th bilayers. We expect that F^ can be 
written as 



F' 



fn ('^n(S„), Cr„, (T„+i) , 



(25) 



where denotes the interaction energy F^ per unit area 
of the n-th bilayer and (i„(S„), defined in Figure El is the 
spacing between the n-th and (n + l)-th bilayers. For 
simplicity we assume that fn depends only on the spac- 
ing dn and the tensions (7„ and an+i H22|l . This implies 
1) = /^(<^n,o-„+i,f7„) since the effect of the 
bilayer curvatures is neglected ^]. Because the bilayer 
cylinders may not be coaxial, we define the inter-bilayer 
spacing (i„(S„) as the shortest distance from the point S„ 
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on the n-th cylinder to the {n 
and tpn defined in Figure 



l)-th cyhnder: With Sn-\ 



1 /2 

sin^ Tpn S„+i 



(r„+i - r„) + 5n+i cosi/i„ 

COS-0n Sin'^ Ipn l^n+l ^ ^ 



'5:^+1 



(26) 



The potential gives rise to an inter-bilayer pressure 



Pn 



dd„ 



(27) 



This pressure may not be equal to the real inter-bilayer 
pressure Vn owing to the fact that the definitions of 
and dn are artificial. Nonetheless, Pn should be a good 
approximation to T'n, especially when \5n+i\ ^ (rn+i— r„). 
For typical myelins, we expect p„ > |16| and moreover, 
Pn should satisfy the force balance equation for the n-th 
bilayer (Laplace's formula) |29| : 



dpi 



-Pn-l - Pn 



Knl + K.n2 



(28) 



where dp^ is the water pressure difference across the bi- 
layer, and Knl and Kn2 are the two principal curvatures 
of the bilayer (see (d)) and Fig. [Tjl. 

4.3 Modeling Myelin Coiling 

The myelin structure defined in this work has three im- 
portant properties: (a) it is postulated to be in quasi- 
equilibrium owing to its slow growth rate (Sect. ^16, 17 , 
(b) the bilayer tension may vary from bilayer to bilayer, 
and (c) the constituent bilayer cylinders can be non- 
coaxial. Given these properties, we will demonstrate that 
coiling may lead to a decrease in the elastic energy F'^ 
()2;-i|l . Therefore, a straight myelin is unstable to coiling if 
the decrease in F'^ is large enough to compensate for the 
net increase in other energies (see ((T7|) . (pn|) . and (j^ l. In 
Section[Slwe will derive the energy changes associated with 
myelin coiling. By analyzing a two-bilayer case in detail in 
Sectional we will show that coiling may indeed lower the 
myelin energy. 

We will employ the two-step deformation developed in 
Section 121 to investigate the energy changes due to myelin 
coiling (Figs. |21and|Hll. For the sake of simplicity, two con- 
straints are imposed. First, the myelin length Lm is fixed 
during the deformation (Fig. 0)). As in the two-dimensional 
case (cf. the text below (QJ), this artificial constraint would 
not invalidate our model. Secondly, we do not allow ex- 
change of material, surfactant or water, between bilayers. 
Because water and bilayers are virtually incompressible un- 
der moderate pressure this constraint implies that the 




Figure 8: Two-step deformation of a myelin. As in Figure 
13 this procedure transforms a myelin from state (i) into 
state (ii) and then into state (iii), with the constraints 
of fixed myelin length and fixed volume for each bilayer 
cylinder. For clarity we only show two adjacent bilayer 
cylinders (Fig. [SJ. Initial state (i): The myelin consists of 
straight, coaxial bilayer cylinders with radii r„ = fn and 
myelin length Lm = Li ~ Li — Ln (Fig. 0J. Intermediate 
state (ii) : The myelin curves while its constituent cylinders 
are still coaxial. In this state f?„ = The constraints 
yield that Rnd = Li, r„ — fn, and therefore An = An 
(see H14f) and \ib\ ). Final state (iii): The cylinders are 
no longer coaxial. The eccentricity vector E„_|_i is defined 
in Figure All the cylinders except for the innermost 
one {n — \) may move toward the center O of bending 
curvature, i.e., 5n+i < 0. In state (iii), ri = fi, Li = Zi, 
Rn+i = Rn + 5n+i < Rn, and therefore L„+i = i?„+i6'_< 
Ln — RnO- The constraints lead to r„ > f„ and An < An 
forn > 1 (see ^ and ^). 
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volume of each bilayer cylinder is conserved during the de- 
formation. (See the Discussion section.) 

Figure IHl gives a detailed description of the two-step de- 
formation. The myelin is deformed from state (i) into state 
(ii) and then into state (iii). In state (i) the myelin is 
straight and its constituent bilayer cylinders are concen- 
tric, i.e., E„ ~ for all n. From (i) to (ii), the myelin 
curves with the condition that the bilayer cylinders are still 
concentric. As a result, the curvature energy F'^ increases 
while the elastic energy and the inter-bilayer interac- 
tion remain unchanged. From (ii) to (iii), the cylinders 
are displaced so that they are no longer coaxial (Fig. [SJ. 
In this step the displacement 8n may become negative for 
n > 1, reducing the bilayer area An and thus the elastic 
energy F° . 

Hereinafter, symbols like fn, Ln, and p„ with an over- 
head bar refer to the initial state (i), where the myelin is 
straight, and those without refer to the final state (iii), 
where the myelin is curved (Fig. (SI. The terms "coil- 
ing" and "bending" are used interchangeably because they 
are energetically equivalent for myelins composed of fluid 
bilayers. Since the myelins considered here have large as- 
pect ratios, their end caps are neglected. However, the end 
caps can be neglected only in terms of energy calculation. 
Their presence is necessary for volume conservation and 
force cancellation (cf. text below (^). 



1/2 



(31) 



where (5i = by construction and 5n < for n > 1. 
Equation (|30(l implies that not only r„+i > f„+i but also 
(r„+i - r„) > (f„+i - f„), i.e., both r„+i and (r„+i - r„) 
increase when the myelin is transformed from state (i) to 
state (iii). The increase in (r„+i — r„) hints the possibility 
of a decrease in the inter-bilayer interaction F^ (|25|l . Merg- 
ing (|25|l and (|30|l leads to the change in the inter-bilayer 
spacing 



where 
Sdni 



Sdn = d„[(iii)] - d„[(i)] 

= d„(S„) - (^'n-l-l - ^^n 
= 5dnl+6dn2, 



(32) 



2 Ri 



n+l ~ ' n) 




3rn I Sj 



5 Bending of a Myelin 

Now we will use the model defined in the previous section 
to investigate myelin coiling. As in Section |31 we will first 
derive the geometrical changes of a myelin undergoing the 
two-step deformation (Fig. |SJ). Then we will determine 
the respective energy changes. For a myelin composed of 
N bilayers, state (iii) is described by N parameters, 
62, . . . , 5n-i, and 5n {Si = 0; see Fig. EJ. 

The two-step deformation is constrained by the condi- 
tions of fixed myelin length Lm = Li = Li = Ln and vol- 
ume conservation Vn ~ Vn for each bilayer cylinder (Sect. 
14.3(1 . With these two constraints and (Figs. EIEI and|HI) 



Rn — Ri 



the length L„, radius r„, and area An of the n-th bilayer 
cylinder in state (iii) are given by (C^ ^'^'^ 03) 



and 



5dr, 



■Jn+1 



COS Ipn 



_ sin^ Ipn Sn+1 
2 fn 
iipn sin^ Ipn (5^+1 



Equation H31() gives rise to the area change 



6An 



AJ(iii)]-A„[(i)] 



An \ ^ Sj An I \ ^ 5j 

~Y^^r'i~ T 1 ^ i?r 



o 



Rf 



< 0, 



(33) 



where = and thus SAi — 0. Equations 
allow us to calculate the energy changes due to the two- 
step deformation, as shown below. 

The curvature energy change (cf. (|20|) ') 



n c 

RnO^Lnll + Y,^ \ < L„ 



Ri 




-1/2 



(29) 
(30) 



SF"" 



F^[(iii)]-F^[(i)] 



N T 

n=l 



2nfnRi 



> 0, 



3 = i 



(34) 
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where L„ = A„/27rr„ = Li and (cf. (|21|l l 



the inter-bilayer pressure of state (i) (cf. (|27|l 'l 



1 H 

^ 4i?? 



(35) 



Pn = 



(38) 



The first-order terms (^ 3nX)^j) H34|) result from the 
second step of the deformation. The exact functional form 
of gn is unimportant because the lowest-order terms (~ 
K„) dominate. This is also intuitively reasonable in that 
when Ri < oo, the major contribution to SF'^ 

should come from the deformation from (i) to (ii) (Fig. IHJ. 
Therefore we will not derive gn here. 

Eauationl33l shows that for n > 1, the area An decreases 
in the two-step deformation. This in turn reduces the elas- 
tic energy F° (see (j221l and (jSHJ)- With the change in 
f is given by 

N 

SF^ EE F^[{ui)]-F^m = Y,SF: 

n=l 

N r 1 



and the spacing change Sdn is given by Ip^ . Using 
and 



S<T„ 



dOn 



bAn 



5Al 



\dAl); 2 



(39) 



we determine the change in the inter-bilayer interaction 
energy as well as the net change in the myelin energy ( H17|l , 
(123, and 



6F' 



n+l 



k=2 j=l 
+ C„5 5n+l 



E 

n=l 
N 



dAn 



2\dAl,, 



5Ai 



O 



(40) 



and 



EAnO'n Sj 

n=l j=l 



6F 



F[(iii)] - F[{i)] = + SF'' + 6F' 



Ar. 



N J 



dOn 

dAr. 



An ^ (Jn 



N 

E 

n=l 
N-1 



J- 

Ri 



O 



EAn 
~9~ 



Sji + ^ Cnk ^ 

k=0 j=l 

3 n+l J ^ (5 - 



(36) 



fc=2 j=l 



where <5Ff — since (Jyli = 0. 

The change in the inter-bilayer interaction F' (|24l) takes 
a complex form because (5yl„ ^ and (5cr„ 7^ for n > 1: 



+ C„5 ^„_|_i 



(41) 



where 



6F' EE ^^'[(iii)]-F'[(i)]= 5]<5F; 

n=l 

= E / C'^n /^(rfn(S„),Cr„,Cr„+i), 



27rf„i?2 2i?2 

(^n+1 - rn)Pn, 



where is defined in l|25|) . Expanding the above equation 
around state (i) results in 



E 

n=l 



-PnSd„ 



dpn\ Sd^ 

ddnJi 2 



TT-^ 5(T„ + — S(Jn+l 



N-1 



E 



(37) 



Bn 

Cnl 
Cn2 

Cn3 



d{AAn 
dAn 

r 1 ( 9pn 



> 0, 



(42) 



/ dcTn- 




\dAn- 


fl/ 


dcTn \ 


An 


dAn); 


\-{ 


OF, 




OA, 



A. 



n+l, 



(43) 
(44) 



n=l 



where the area element (Fig. [Tj) 



dAn = RldO ■ Tndlpn \ "'^ + ^ ^" 



5i = by definition, and ^„ < for n > \. In spite of the 
fact that equation ||1T|| allows 5„ to be positive or negative, 
we are only interested in the cases where (5„ < (Figs. |S1 
andlH)). The terms Bn and C„o come from the curvature 
energy F'^ (f^ . C„i from the elastic energy F° ((23), and 
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Cn2 — Cn5 from the inter-bilayer interaction energy F' H24|) . 
As mentioned in the text below (|^ . C„o is neghgible. In 
addition, we expect p„ > and (dpn/ddn) < for myehns 
[T?)j . Equation H4HI seems too compUcated to be useful 
for investigating myelin coiling: there are 2N initial condi- 
tions, an and f„, and N independent parameters, i?i and 
5n {Si = 0; see Fig. 0). However, this equation is not as 
worthless as it seems. In order for a myelin to coil, SF 
(|41|l has to be less than zero. The two-dimensional case 
of Section O suggests that if some of the bilayer tensions 
CT„ are sufficiently large, the decrease in F° might 
be large enough to yield a negative SF. In the next sec- 
tion we will show that the proposed mechanism of myelin 
coiling can be easily understood with the special case of 
N — 2. The enormous number of parameters in equation 
(|41|l merely reflects the complexity of myelinic structures: 
for instance, two apparently similar myelins may have very 
different bilayer tension profiles and therefore one may coil 
and the other may not. 



6 Bending of a Two-bilayer Myelin 



In this section we inspect closely a special case of myelin 
coiling, where the myelin is composed of only two bilay- 
ers, i.e., N = 2. One may think of this two-bilayer 
myelin as a generalization of the two-dimensional tube of 
Section O Besides, equation H41I) suggests that in terms 
of energy calculation, an A^-bilayer myelin can be viewed 
as a multilamellar tube composed of {N — 1) nested two- 
bilayer myelins. Therefore the results of this two-bilayer 
case should be applicable to multiple-bilayer myelins. In 
the following we will first derive a simple criterion of coiling 
instability for a two-bilayer myelin and study a numerical 
example in detail in Sections 16.11 and 16.21 Then in Sec- 
tion |^31wc will show how this two-bilayer example can be 
generalized for myelins of many bilayers 

For a two-bilayer myelin, two independent parameters 
Ri and S2 determine state (iii) (Figs. [Sland|Sl). Equation 
(|41|l can be rewritten as 



SF ^ i^[(iii)] - ^^[(i)] 



2 



B+-^J2^k+C5S 



Ri 



k=0 



(45) 



where ^2 < (Figs. |31and|Hl), 



B 

Co 



B1 + B2 



A. 



Kl + K2 



Ai 2TTflR{ 



>o, 



A2 



= 92- 



r2 
) . 



= C21^ = CT2— > 0, 

Ai n 
= C12 = (r2 - ri)pi > 0, 



C2 
C3 



(46) 
(47) 

(48) 

(49) 



and Ci4 H43|l does not exist because SAi = H33|l and Sai = 
(1(23 and (ISni)- Although equation (gSJ does not dictate 
the sign of (52, we only consider the case of ^2 < because of 
the two-step deformation (Figs. |Sland|SI). The terms B and 
Co originate from the curvature energy F"^ (I20() , Ci from the 
elastic energy F° (I23|l . and C2, C3, and C5 from the inter- 
bilayer interaction energy F^ (124(1 . Here we assume that the 
outer bilayer is under tension ((T2 > 0) and the inter-bilayer 
interaction is repulsive {pi > 0). We expect C3 < and 
C5 > because the inter-bilayer repulsion increases as the 
bilayer tension a2 or the inter-bilayer spacing di decreases, 
i.e., df[/da2 < and {dpi/ddi) < 0. Therefore, Ci and C2 
contribute to coiling, while C3 and C5 hinder coiling. 

6.1 Criterion of Coiling Instability 

Equation (|45|l gives rise to a simple stability criterion when 
all the 0((5|) terms are negligible: a straight, two-bilayer 
myelin is unstable against bending or coiling if there exists 
a set of {S2, Ri) such that 



SFr^B- 



Ri 



[Ci + C2 + C3] < 0, 



(50) 



where Co is dropped because B ^ \CoS2/Ri\ (see the text 
below H35(l 'l. By setting SF = we can define a threshold 
of S2, denoted as ^th, for coiling instability. If one can 
find a finite i?i and a negative S2 of magnitude larger than 
I (5th I, then state (iii) is energetically more favorable than 
state (i) (Fig. [HJ. However, S2 and i?i must be bounded 
by the geometry of the myelin: for example, 1^2 1 must be 
less than (f2 — ^i), and Ri must be greater than r2. The 
equilibrium value of (^2,^1) can be determined when the 
terms of C((5|) in 145(1 are known. 

If the tension (T2 is large enough so that Ci » IC2 + C3I, 
i.e., the inter-bilayer interaction is negligible, criterion ((50(1 
is reduced to 

SF^B+^Ci<0. (51) 
Ri 

Furthermore, when f2 ~ f 1 ^ (f2 ~ '^i); we expect B « 
k/RI (see (|S5)l and (03)) and Ci « (T2. In such case the i52 
threshold for coiling is given by 



<5th = 



CT2i?l 



(52) 



Now we investigate a two-bilayer myelin composed of typ- 
ical lipid bilayers with parameters 



r-2 



20 /xm (Fig. HJ, 
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Ri 

K 
Pi 



f2-fi= 2.4 nm gT], 
30 Aim IHlinillSI, 



= 10^ 



erg [11, 



dyne/cm^ (Eq. 
lO^dync/cm^ (Eq. ^) 



and the temperature T = 4.1 x 10^^^ erg (300 K). For sim- 
plicity the water pressure difference dp^ is assumed to be 
zero. The values of these parameters are carefully chosen 
such that the results of this example can be readily gener- 
alized for myelins of many bilayers (Sect. I6.3|l . Equation 
(|28|l with the above parameters implies Ci ~ 02 ~ r2pi — 2 
dyne/cm J^, Jfi, and thus Sth ~ —0.02 A according to H52|l . 
Because \6th\ is much less than typical bilayer thickness or 
inter-bilayer spacing (Fig. ^ 57^, there exists a negative 
S2 such that I 1 > \Sth\- Therefore state (iii) is more sta- 
ble than state (i) (see (I50|l and Fig. (HJ. This numerical 
example suggests that criterion (|51|l can be met easily and 
a two-bilayer myelin may bend or coil if the lateral ten- 
sion of the outer bilayer is sufficiently large. However, two 
questions arise with regard to the above example. Can we 
ignore the 0(62) terms in equation H45(l ? Moreover, is the 
inter-bilayer interaction negligible, i.e., is it possible that 
Ci ^ |C2-f C3I? Because the complete functional form of /,\ 
(|25|l is not known, it seems difficult, if not impossible, to 
answer these two questions with logical rigor. In the next 
subsection we will try to address these issues by means of 
some physical arguments. 

In fact, equation (|52|l can also be viewed from another 
perspective: it provides a lower bound for (T2 in a coiled 
myelin when |(Sth| is replaced with the inter-bilayer spacing, 
i.e., a coiled two-bilayer myelin with coiling curvature l/i?i 
(|13|l must have a (T2 greater than K/{diRi). The actual 
5'2 is expected to be larger than this value because of the 
contributions from the inter-bilayer interaction, the higher- 
order terms in Ef^ l(T^ . and the 0(62) terms in In 
principle, we can obtain a more sophisticated lower bound 
for (T2 using equation (|45|l if the functional form of /{ is 
known. 

6.2 Validity of the Criterion 

Equation (|52(l is valid only when (a) f2 « fi 3> (^2 — ^1), 
(b) Ci ^ IC2 -I-C3I, i.e., the inter-bilayer interaction is neg- 
ligible, and (c) the 0(i5|) terms of can be ignored. 
Condition (a) is obviously met by the numerical example 
given below (|52l) . Intuitively, condition (b) would be sat- 
isfied if the bilayer tension CT2 is sufficiently large. For a 
straight two-bilayer myelin, a2 is balanced by the net pres- 
sure difference Ap across the outer bilayer cylinder H28|l : 
(72 = r2^P ~ ^2{dp^ where pi (|38|l represents the 

strength of the inter-bilayer interaction. This force balance 
equation implies that a2 increases with f2 or Ap. Because 
Pi (> 0) tends to decrease as a2 increases, we expect a 
sufficiently large f2 or dp2 would yield a CT2 large enough 



so that the inter-bilayer interaction can be neglected. In 
the following we will examine closely whether conditions 
(b) and (c) are satisfied by the example below H52|l . Our 
strategy is to compare the orders of magnitude of the terms 
of interest. 

Equations ^ and with 5th ~ -O.O2A yield 



'th| 



6di 



SA2 
A2 



(53) 



Here the use of symbol means we are only interested 
in the orders of magnitude of the terms on its both sides. 
The above equation in turn implies that the 0(i5f ) terms 
in 6di (|32|l are negligible. Experimental results suggest 
that the bilayer tension varies linearly with the fractional 
change in bilayer area when it is sufficiently large (> 0.5 
dyne/cm for typical lipid bilayers) |25[ I26| . Thus, in our 
case 



dA2Ji- 



Ka, 



(54) 



where K^. is the area compressibility modulus of the bi- 
layer. For typical lipid bilayer membranes, if a ~ 200 
dyne/cm [221 ■ Equations H39I) and (|54|) imply that the 
fractional change in (T2 is also small: 



Sa2 

CT2 



KaSA2 
02 A2 



(55) 



where (T2 = 2 dyne/cm. By virtue of H53|) and (|54|l . it is 
straightforward to show that in JF^ (j^Hl, the 0{Sj/Rl) 
term is much less than the 0(62/ Ri) term and thus can be 
ignored. 

The change in the inter-bilayer interaction potential /{ 
reads (see and ((23) 



fc=2 



(56) 



where, with and 



V2 = 
V-i = 

= 

= 

-De = 



-Pi5di, 



5al 



2-/i 2 



(dpA Sd]_ 
\ddJi 2 ' 



and Sdi ~ |i5th| ~ 0.02 A in our case. The first-order terms 
V2 and 2?3 lead to C2 l|T7|) and C3 I^HJl, respectively (see 
^7^). The second-order terms V^, V5, and Vq contribute 
to the 0(62) terms in H45|l . In our case Ci « 2 erg/cm^ 
H46|l is certainly much larger than C2 = 2.4 x 10~* erg/cm^ 
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(|77jl because » (r2 — ri) and (T2 ~ ''2j'i (EH}- Therefore Given "K 
C2 and I?2 can be neglected. Besides, we expect that 



10 ^, (72 w 2 dyne/cm, if a ~ 200 dyne/cm, 
1 nm, equations (|^ . (j^SJl, 



/i~PiA and 



9pi 
ddi 



Pi 
A ^ 



(57) 



where the length A can be as small as the decay length 
(2 - 3 A) of hydration pressure |2Z1 or as large as the inter- 
bilayer spacing (several nanometers) if the Helfrich repul- 
sion j28j dominates. This means in our case 



Sdi 



\s. 



th 



A 



10" 



10"^ < 1 



(58) 



and |(5th|Ci/i?i 



C5(5jjj ~ 10 ^ dyne/cm (see 



lO^'' dyne/cm is much larger than ~ 
3))). Hence we can ignore 
those 0(62) terms of ()45|) arising from D5 . The functional 
dependence of /{ on (T2 is necessary in order to decide 
whether P3, 2?4, and Vq are negligible. To go further, 
we consider two extreme situations. The actual /i(o'2) is 
expected to fall in between these two extremes. 

In the first extreme situation, we assume 172 is large 
enough so that the bilayer thermal undulations P[2| can 
be ignored. This assumption implies that tT2 is coupled 
with fl mainly via surfactant number density pg (number 
of surfactant molecules per unit bilayer area) of the outer 
bilayer. Intuitively, when (T2 is sufficiently large (close to 
the lysis tension) j^, /{ should be roughly proportional to 
Ps- With fixed inter-bilayer spacing di (see (|56|) '). equations 
(|^ . (|^ . and the assumption of /{ oc ps lead to 



Spi 5J[ 
Pi ^ fi 



SA2 ^ (72 6(72 
M A'a (72 



and thus 



where 



9(72 602 (72 ' 



(59) 



(60) 



y _ Sf[/fl Sf{/f{ (72 _ (72 ^ ^g^^ 

5(72/(72 SA2/A2KA Ka 

In our case a2 = 2 dyne/cm, Ka « 200 dync/cm, and 
hence Y ~ 10"^. By use of H57|l . H59() . and (|60|1 . we can 
rewrite 2?3, 2?4, and T^q: 



(72 \ (72 



Y(Y - l)piA 



-SpiSdi 



5(72 
0-2 

6pi \ ( Sdi 



Pi 



A 



pi = 10 dync/cm , and A 
and (|58|l imply 

\SF''\ _ (72\6A2\ al \5a2\ 



A2 



A2 



Ka (72 



> IP3I > \'Di\ or \Vq 



Therefore, C3 H48|l . which arises from 1)3, as well as those 
0{52) terms of (|45f) contributed by and are negli- 
gible in the example below H52() . We also note that with 

m, in, nzi), (Eoi, and inu. 



C3 



9(72 

Pi A 



fi 

Y^Ka 



Ka 

i CT2 

- 10"^ - 10""* 



-YKa 



CT2 



erg/cm^ 



which has approximately the same order of magnitude as 
C2 = 2.4 X 10^^ erg/cm^ (I47|l . This means in our case 
coiling might even lead to (C2 -I-C3) > 0, i.e., a decrease in 
the inter-bilayer interaction (see (|^ . (1^ . (|5n|) . and 
the text below (jl^ ). 

Now we consider the other extreme situation: the de- 
pendence of /{ on (72 originates from bilayer thermal un- 
dulations |19l 1301 m] . Here we adopt the functional form 
derived by Seifert (Eq. (9) of [^01) and express /} as 



/i = /i(di) + 



y 



ndl sinh^(i/)' 



(62) 



where b = 0.1, di = [7=2 - fi), and y = {a2/hTY/^di/2. 
The term f\ is the part of /{ that is independent of (72- 
The last term is just the Helfrich potential with variable 
membrane tension. The dependence on ai is neglected 
because cti is unchanged in the two-step deformation. With 
H62(l and the parameters given below (|52|l , straightforward 
calculations show that Ci w 2 erg/cm^ H46(l is much larger 
than IC3I « 0.04 erg/cm^ ^ and that both P4 and 
are much less than V^. Therefore, C3 as well as the ©((51) 
terms of associated with I?4 and are negligible. 

In the previous subsection we use equation (|52(l to in- 
vestigate the coiling instability of a two-bilayer myelin (see 
(|51|) '). This equation is valid only when (a) f2 ~ ri ^ 
(^2 — ^1), (b) Ci » IC2 -I-C3I, i.e., the inter-bilayer interac- 
tion is negligible, and (c) the ©((^j) terms in can be 
ignored. We are interested in the example described below 
1)52(1 . which clearly satisfies condition (a). In this subsec- 
tion we examine conditions (b) and (c) by estimating the 
orders of magnitude of the relevant terms. We show that 
these two conditions are also satisfied. Thus, equation ((5^ 
is indeed applicable to the case of interest. 

6.3 Generalization for Multiple-bilayer 
Myelins 

In principle the results of Section 16.11 can be applied to 
myelins of multiple bilayers. This is because in terms of 
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energy (see (|17|l and H41|l 'l. an 7V-bilayer myelin can be 
decomposed into {N — 1) nested two-bilayer myelins. Al- 
though the bilayer tension f7„ ()22(l may vary from bilayer 
to bilayer, equation H41|) implies that the coiled state is still 
favorable as long as the net decrease in F° is larger than 
the net increase in + F\ 

Equation 141(1 and its implications for myelin coiling can 
also be understood with a gedanken experiment, where 
a multiple-bilayer myelin is constructed from two-bilayer 
myelins. We take the numerical example below ()52|l as our 
initial two-bilayer myelin. The bending curvature 1/Ri 
and the eccentricity 62 are allowed to vary freely so that 
the myelin is in equilibrium. The equilibrium state has a fi- 
nite Ri and a ^2 of magnitude larger than \Sth\ i~ 0.02A). 
This means 6F (|45|l is less than zero. Put it in another 
way, the two-step deformation (Fig. ISj) produces a surplus 
of energy. Since \6th\ of this initial myelin is quite small, 
the equilibrium value of \d2\ is expected to be much larger 
than \Sth\- Therefore the energy surplus \6F\ should also 
be much larger than the curvature energy cost of coiling 
SF'^ ((34|l . We keep i?i unchanged and insert another two- 
bilayer myelin into the inside of the initial one. Now the 
system becomes a four-bilayer myelin. The energy sur- 
plus \SF\ may increase or decrease as a result of this new 
addition. We can repeat this process and add more two- 
bilayer myelins to the system as long as the energy sur- 
plus is not depleted, i.e., 6F is still less than zero. Even 
when 6F > 0, we can increase the bilayer tensions in those 
existing two-bilayer myelins so that 6F becomes negative 
again. The above gedanken experiment shows clearly that 
a multi-bilayer myelin can coil to lower its energy when the 
bilayer tension is sufficiently large. 

Let us consider the case depicted by Figure Eb- In this 
case the material influx density for the core is increasingly 
greater than that for the shell. Eventually the bilayer ten- 
sion in the shell will be large enough to cause the myelin 
to bend or coil. To be more specific, the free energy anal- 
ysis ((4111 can be done in the following manner. First, the 
eccentricities (5„ (Fig. EJ of the bilayers in the core are ar- 
tificially constrained to be zero, while those in the shell are 
free to vary. We start with state (i) (Fig. |SJ), i.e., i?i — > 00 
and 6n = for all n. By symmetry state (i) is the equi- 
librium configuration for straight myelins. Secondly, we 
perform the two-step deformation (Fig. ^ and allow the 
system to reach equilibrium at finite i?i . This means those 
variable eccentricities 5„ are at the values minimizing the 
total energy F lfT7|l. Because every (5„ in the core is set 
to zero, the deformation of the core only costs curvature 
energy ((34(1 . Previous analysis of the two-bilayer myelin 
suggests that when the bilayer tensions cr„ in the shell are 
large enough, the two-step deformation would lead to a net 
decrease in F and thus state (iii) becomes energetically fa- 
vorable. Lastly, the constraint of (5„ = in the core is 
dismissed and the system is allowed to equilibrate again. 
Obviously, if state (iii) has been more stable than state (i) 



as a result of large bilayer tension in the shell, it would be 
still favorable because the new degrees of freedom would 
only lower the total energy F further. Therefore, we con- 
clude that the myelin of Figure[5j3 would bend or coil when 
the bilayer tension in the shell is sufficiently large. In the 
sense of coarse-graining, this case can be treated as a "two- 
bilayer" myelin, whose inner and outer bilayers correspond 
to the core and the shell, respectively. 

7 Discussion 

Our model predicts that a sufficiently large bilayer tension 
may cause a myelin to coil. The coiling mechanism is ex- 
plained with the two-bilayer case studied in detail in the 
previous section. In the following we will discuss the exper- 
imental techniques that might be used to test our model. 
We will also look into some important features as well as 
implications of our model. 

In order to test the proposed model of myelin coiling, 
one needs to be able to either measure or manipulate the 
bilayer tension of a myelin. Optical tweezers may be a 
promising tool [62i . When a bilayer is optically tweez- 
ered, surfactant molecules are pulled into the trap set by 
the tweezers and thus the bilayer tension increases. An ad- 
vantage of optical tweezers is that they can penetrate into a 
myelin and grab multiple bilayers. However, laser-induced 
tension is usually weak (10"'^ — 10~^ dyne/cm) [311. K 
remains a question whether such a weak tension can in- 
duce myelin coiling. Another method is micropipet suction 
(251 ■ In contrast to optical tweezers, this method can 
produce a bilayer tension as large as the lysis tension . 
Nonetheless, it can be applied only from the outside of a 
myelin. Although suction of a myelin with a micropipet 
is an interesting problem in its own right, whether it can 
cause myelin coiling is still uncertain. 

The second constraint described in Section 14.31 requires 
volume conservation for each bilayer cylinder. It should be 
a fair condition for defect-free myelins j35( since water per- 
meabilities of typical lipid bilayers are quite small ^1 EEl ■ 
This constraint is also necessary in our model because it 
eliminates a potential inconsistency between properties (a) 
and (b) described in Section IH^ if bilayers can exchange 
material, a nonuniform bilayer tension may induce a suf- 
ficiently large material flow between bilayers to relax the 
tension that drives coiling. 

For simplicity the spontaneous curvature Kq ((18() is set 
to zero in this work. A nonzero or nonuniform Kq would 
only alter the curvature energy F"^ ((20(1 by some finite 
amount and therefore does not weaken the proposed coiling 
mechanism at all ^JE]- In addition, the cross sections 
of the bilayer cylinders are assumed to be circular (Fig. 
O. Although this might be a fair assumption for straight 
myelins, there is no reason why this is still the case for 
coiled myelins. Because (5„ is in general nonzero in state 
(iii) (Fig. IH)), the shape of the cross section should change 
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accordingly. However, allowing this to happen would only 
lower the energy of a coiled myelin further (cf. H41|) ') and 
thus would not render our model invalid. 



The myelin free energy F H17|) is the sum of the curva- 
ture energy F'^ , the elastic energy F° of bilayer tension, 
and the inter- bilayer interaction energy F^. Because F'^ is 
known (|2U|I . determining the change in F'^ due to coiling 
is straightforward H34|l . The change in to lowest order, 
is given by ^(T„(5yl„, where (t„ is taken as an indepen- 
dent parameter (see (|36|l '). Thus, direct knowledge of F° is 
not crucial for studying myelin coiling at the level of this 
work. This, however, is not the case for F^ l|24|l . In order 
to analyze the coiling instability in detail, the functional 
form of (|25|) is required even though we have shown that 
its contributions can be neglected under certain conditions 
(Sect. EJ. Yet, to our knowledge, the complete functional 
expression of /,\ is still not known. Therefore more inves- 
tigations of are needed. In the context of this work, /,\ 
has to satisfy the force balance equation if^ . 



In the example below H52|l we use a small inter-bilayer 
pressure (pi = 10^ dyne/cm^) while our previous work |16| 
suggests that the inter-bilayer pressure (|?7jl in a multi- 
bilayer myelin is, on the average, larger than this value. 
This, however, is not a contradiction. In view of equa- 
tion H28|l . a two-bilayer myelin with (T2 = 2 dyne/cm, 
r2 = 20 /zm, and dp!^ = dyne/cm^ must have an inter- 
bilayer pressure of about 10^ dyne/cm^. As for myelins 
composed of multiple bilayers, equation also allows 
a much larger inter-bilayer pressure p„, provided that the 
net pressure difference {dp^ + pn-i — Pn) yields a reason- 
able bilayer tension (t„. Thus, a multi-bilayer myelin may 
have a large inter-bilayer pressure Pn as long as the net 
pressure gradient across the bilayers is small. This leads 
to an interesting question in the physics of myelins: what 
determines the width (diameter) of myelins? Our model 
suggests bilayer tension may be an important factor. As 
just mentioned, an inter-bilayer pressure of 10^ dyne/cm^ 
(« 10~^ atm) with a radius of 20 fim yields a bilayer ten- 
sion of 2 dyne/cm, which is close to the lysis tension of 
lipid bilayer membranes (1-20 dyne/cm) [2SI- In contact 
experiments myelins of diameter 10 - 50 /Ltm form immedi- 
ately after bulk surfactant is contacted with water (Sect. 
13) [HJITIIS]. The surfactant concentration at the contact in- 
terface should be sufficiently high so that the inter-bilayer 
interaction is repulsive. Therefore, even in view of the 
complex kinetics of surfactant dissolution |H7j , an obvious 
upper bound for myelin diameter is prescribed by the lysis 
tension of the bilayer and the outermost inter-bilayer pres- 
sure in myelins. By the same token, we can estimate the 
outermost inter-bilayer pressure using the observed myelin 
diameter and lysis tension. 



8 Conclusions 

In this work we propose a model to explain coiling of 
myelins often seen in contact experiments. This model has 
two important features: the constituent bilayer cylinders 
of a myelin can be non-coaxial and the bilayer tension can 
vary from bilayer to bilayer. The analysis indicates the 
importance of bilayer tension in determining the morphol- 
ogy of a myelin: a myelin can coil if the bilayer tension is 
large enough. The coiling mechanism is in a sense similar 
to the classical Euler buckling of a thin elastic rod subject 
to axial compression |2()j . 

Our model suggests that myelins have complex internal 
structures. The internal stresses such as bilayer lateral ten- 
sion may be very different for seemingly like myelins. As a 
result, our model should not be deemed inconsistent with 
other models |12lli;il[T7j . On the contrary, the existence of 
all these models, including ours, reflects the complexity of 
myelinic structures and suggests multiple routes to myelin 
coiling. However, the following questions are also legiti- 
mate. Is our model applicable to the systems that were 
previously explained with different models jlOl 1121 IT^ ? If 
so, how can we determine which model better explains the 
experimental results? Moreover, owing to the complex- 
ity of myelins, experimentalists are faced with a serious 
challenge: how can one produce myelins in a controllable 
fashion? To answer these questions, more experimental as 
well as theoretical investigations are needed. 
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